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Abstract 

We study the conductivity of a strongly coupled striped superconductor using gauge/gravity 
duality (holography). The study is done analytically, in the large modulation regime. We show 
that the optical conductivity is inhomogeneous but isotropic at low temperatures. Near but below 
^■i the critical temperature, we calculate the conductivity analytically at small frequency lo and find 



it to be both inhomogeneous and anisotropic. The anisotropy is imaginary and scales like 1/uj. We 
also calculate analytically the speed of the second sound and the thermodynamic susceptibility. 
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I. INTRODUCTION 



Recently, there has been a flurry of activities applying gauge/gravity duality to study 
strongly coupled condensed matter systems (for a review, see for example, Ref. A 
large portion of these activities revolves around understanding the so-called high-T c or high- 
temperature superconductors. The high-T c superconductors, such as cuprates and iron pnic- 
tides, are interesting, because not only do they exhibit strong correlation, but they also 
exhibit the existence of competing orders that are related to the breaking of the lattice sym- 
metries. At first, these orders seem to be unrelated to superconductivity, however, a study 
of the effect of inhomogeneity of the pairing interaction in a weakly coupled BCS system 
[2I as well as numerical studies of Hubbard models (3), 4] suggest that inhomogeneity might 
play a role in high-T c superconductivity. Furthermore, the recent discovery of transport 
anomalies in La2- x Ba z Cu04, which are particularly prominent for x = 1/8 might be 
explained under the assumption that this cuprate is a superconductor with a unidirectional 
charge density wave, i.e., a "striped" superconductor pj. Other studies using mean- field 
theory have also shown that unlike the homogeneous superconductor, the^striped supercon- 
ductor exhibits the existence of a Fermi surface in the ordered phase 



sensitivity to quenched disorder 



BE 



and its complex 



In gauge/gravity duality, the strongly coupled condensed matter systems are mapped to 
a weakly coupled gravitational theory on a spacetime with negative cosmological constant 
(asymptotically anti de-Sitter (AdS) spacetimes). In particular, the study of high-T c super- 
conductors is mapped to a study of Einstein-Maxwell-scalar theory of four-dimensional AdS 
black holes, and the superconducting phase transition is understood as the (unstable) black 



ho 
9 



es forming scalar 'hair.' To study the effect of inhomogeneity at strong coupling, Refs. 



ll| introduced a modulated chemical potential, which in the language of gauge/gravity 



duality is translated into a modulated boundary value for the electrostatic potential in the 



AdS Schwarzschild black hole. In particular, Ref. [10[ analytically studied the critical 
temperature and the behavior of the order parameter (condensate) as a function of the 
modulation wave number. 

In this article, we study the conductivity of the holographic striped superconductor by 



adding perturbations to the solution found in Ref. 



10] . Our main goal is two-fold: 



(a) To look for anisotropy and study its scaling behavior, i.e., dependence on the frequency. 



As this is something that can be measured in experiment, studying the anisotropy of 
conductivity can lead us toward a more realistic model of high temperature supercon- 
ductors using holographic techniques. 

(b) To establish that indeed this model exhibits a superconducting transition. 

Following Ref. [u]], we focus on the large modulation regime for an analytical study. A more 
thorough analysis requires a combination of analytic and numerical techniques and will be 
taken up elsewhere. Our study is conducted at temperatures near the critical temperature 
T Rj T c and at low temperature T <C T c . In both regimes, we find that the dynamical 
conductivity does not only depend on the frequency and wavelength of the electric field 
applied to the strongly coupled system, but also has a spatial dependence. We also find 
that at T < T c , at zero wavelength, the real part of the conductivity contains a term 
that is proportional to 8(w), where u is the frequency of the applied electric field. This 
indeed shows that the system exhibits superconductivity below the critical temperature. 
Furthermore, at T < T c , we find that the optical conductivity is anisotropic. At small 
frequency, the anisotropy is imaginary and scales like At high frequency, a numerical 
calculation is required. We shall report on the high frequency regime in the near future. 
At low temperature, we study the anisotropy at A = 1, where our analytic technique is 
applicable, and find the conductivity to be isotropic. Lastly, at T ~ T c , we also calculate 
the speed of the second sound and the thermodynamic susceptibility. 

Our discussion is organized as follows. In Section [Til we se t up the system and fix our 
notation. In Section [TTTl we obtain analytic results for the conductivity at low temperatures. 
In Section [IV] ,we calculate the conductivity analytically near the critical temperature. It 
should be noted that the analytical techniques used in the two temperature regimes discussed 
in Sections II III and IIV[ respectively, differ from each other. Finally, Section [V] contains our 
concluding remarks. 



II. SET-UP 



For completeness, let us briefly review the set-up and the results of Ref. lOj. The 
minimal requirement to study strongly coupled superconductor using gauge/gravity duality 
is to have a scalar field and a U (1) gauge field living in a spacetime with negative cosmological 



3 



constant. The scalar field is dual to a superconducting order parameter, i.e., the condensate, 
while the U(l) gauge field is dual to the current in the condensed matter system. 

To study the strong coupling regime of the superconductor, we only need to study the 
gravity theory at the classical level. In particular, we are interested in finding solutions to 
the classical equations of motion whose boundary values are related to the parameters of 
the superconductor. Neglecting the backreaction of the scalar field and the gauge field onto 
the metric, the problem of studying a strongly coupled superconductor with an effective 
two spatial dimensions is reduced to studying the dynamics of a scalar field with mass m, 
coupled to a U{1) gauge field in the background of a 3+1-dimensional AdS Schwarzschild 
black hole with the following metric 

^ = ^(_ ftW ^ + _^!_ +(i /), a) 

with 

h = 1 - z 3 , (2) 

in units in which the AdS radius is unity. In this coordinate system, the boundary is at 
z = 0, while the horizon is at z = 1. 

We introduce inhomogeneity by including a modulated chemical potential, or in other 
words, the charge density wave in the system is being sourced by a modulated chemical 
potential. The electrostatic potential is then given by 

A t (z,x) = fi(l-5)A < t 0) {z) + fi5A ( t 1) cosQx. (3) 

We are particularly interested in the regime where the modulation wave number Q is a lot 
larger than the other scales in the system, i.e., T c and /i. In such regime, Ap is exponentially 
suppressed in the bulk. 

At T > T c , the order parameter, and thus the scalar field, vanishes, but below T c , we 
have non- vanishing scalar field 

^ = (0 A ) ° ^ yv F(n) ^ ^ 
V 2 n>0 

with F(°'(z = 0) = 1, such that the order parameter is 

(O a ) = J2(°A) [n) cosnQx, (5) 

n>0 



with 

(0 A > (n) = (e> A ) (0) F (n) (z = 0). (6) 
Here, the gap [(O^)] 1 ^ is measured in the unit of the temperature. Furthermore, A = 



A± = | ± y | + m? is the dimension of the superconducting order parameter Oa, with m 
being the mass of the scalar field dual to the order parameter. We shall examine the range 

- < A < -, (7) 
2 2' v ; 

for A = A_, corresponding to masses in the range > m 2 > — |. At large Q, the higher 

modes of the scalar order parameter is suppressed compared to the zero mode 

(o A )(°) - w 2[n/2] / ' 

where [n] denotes the smallest integer > n. Furthermore, when the homogeneous part of 
the chemical potential vanishes , i.e., 5 = 1, the odd modes vanish. 

We can now consider adding the vector potential as a perturbation on top of this back- 
ground. From the on-shell Maxwell action, we can then read the current-current retarded 
correlation which will then give us the conductivity. 

III. LOW TEMPERATURE 

A. Conductivity perpendicular to the direction of the stripes 

First, let us study the conductivity a y , which is perpendicular to the direction of the 
stripes, at low temperature. To do so, we introduce a small perturbation 

\ — j^0a y {z^k)e^. (9) 

The leading order of the equation of motion is then given by 

d z (hd z a y ) + ^ay = (V + (O A ) {0)2 s^FW 8 ) a y , (10) 

with an incoming!! boundary condition at the horizon. We can rewrite this in a Schrodinger 



form 

)2 



d Zt a y + Va y = uj a y , (11) 



1 "Incoming" here means that any perturbation introduced in the spacetime will fall into the horizon of the 
black hole. This boundary condition corresponds to the retarded correlators of the boundary theory. 
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where the "potential" is given by 

V = V (0) + k 2 V 1 , (12) 

with 

V® = (0 A )^ 2 z 2 ^hF^ 2 , Vi = h, 
and we have introduced the tortoise coordinate 

f z dz 1 , (1 - z 3 ) If ,l + 2,z tt\ 

Near the boundary, z* « z — )■ 0, and the horizon is at — > oo. 

We can solve the equation of motion by replacing the potential V with its average (V) 
in a self-consistent manner. The solution with the correct boundary condition is then given 
by 

a y a eV"M<V* 2 (14) 

Substituting the solution into the on-shell action, we obtain the retarded current- current 
correlator 

( J y (u, k) J y (-cu, -k)) = i\Ju 2 - (V^ ] ) -k 2 (V^, (15) 
from which, by using the Kubo formula, we obtain the conductivity, 



(16) 



Here, 



_ Jo dz* Vq 0) \a y \ 



POO , I |2 

Jo z * \ a v\ 

T(2A - 1) (O a )(°) 2 Ui^^-iV^-k^V^ ^ , (17) 



where the integrals are evaluated by giving u an imaginary part, which we will set to zero at 
the end, and by approximating the integrals by evaluating the integrand near the boundary 
We also have 

W ^ ^l, (18) 
Jo z * \ a y\ 

and thus 

9 / I \ 2 ( 1 - A ) 



(C0 = r(2A-l)(O A )^ [2i^u 2 -(Vr)-k 2 \ . (19) 



The approximation we used to evaluate the integrals renders this approach insensitive 
to the details of the horizon in the bulk, which corresponds to the infrared regime in the 
boundary theory. Therefore, this approach is valid when the gap is a lot larger than the tem- 
perature, [(Oa)^} 1 ^ S> 1 (in our notation, the gap is measured in the unit of temperature, 
see Section HB . 



1. A = 1 



Eq. [19] can be solved exactly for A = 1 and we obtain 



a v (u) 



(0)' 



(20) 



Now, let us include the first order correction. For 8 ^ 1, we have 



hd z a y (u, fc)] + ^a y (oo, k) = [k 2 + (O A ) (0)2 2 2(A " 1) F (0)2 ) a y (u, k) 



+(O A ) (0) ^-DF^fW [ay(u>, k-Q) + a y (u, k + Q)} 

+ (O A ) (0)2 Z 2(A-1) F (0) F (2) fc _ 2g) + ^ fe + 2 g ) ] 

(21) 



Since and F 1 ^ have Q-dependence suppression, we can solve this perturbatively 

a y (z; u, k) = a y (z; u, k) + a y (z; u, k), 



(22) 



with boundary condition 

a y (z = 0; to, k) = a y (z = 0; u, k). 
For A = 1 this yields at zeroth order 



al(z; to, k) = a° y (z = 0; to, k) e iQkZ * , tt k = Jtu 2 - (OJW - k 2 



(23) 



(24) 



and at first order, 



(C*i) (0)2 , (i),a°Jz = 0;u,k + Q) . . . 

\ 11 (Vp ) — (g^fc + Q 2 . _ g^feZ*^ 



Q + 2k 



{O 1 )^ 2 nA2) y y (z = 0;u J ,k + 2Q) 



4 g W> 
-(Q -+ -Q) , 







■) 



(25) 
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where 



(V rW > = JZ°dz.F®{l-3?) \a y \ 



(Oi)(P) 



i = l,2 , 



(26) 



and we have approximated the integral by evaluating the integrand near the boundary. 
The on-shell action then becomes 



doj dk dk' 
(2tt) 3 



dy(z = 0; to, k)iQ(co; k, k')a y (z = 0; —u, k') 



(27) 



where 



g(u; k, k') = n k 5(k> + k)+ {<Dl) ^ 0l) gUfa k, k') + W ^ 9 {2 \^ k, k') , (28) 
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nk ~ Q " + k-Q) + Qk t Q ~? k Kk' + k + Q) 



Q-2k K ^' Q + 2k 



+ A; + 2Q) 



Q - fc Q + k 

from which we deduce the retarded current-current correlation function 

(J(-w,x) J(a;,fc)> = — / —e- ik *g(uo;k,k') . 

Ztx I In 



(29) 



(30) 



The dynamical conductivity is given by the "pre- averaging" Kubo formula (see, e.g., 12]) 



-ikx 



a y (uj, k, x) = (2tt) — — (J(—oj, x) J (to, k)) . 



IUJ 



(31) 



Explicitly, 



a y (tJ, k, x) 



1 



(0)//-O \(2) 



Q 



4Q 



(32) 



^fc-Q - jQ x ^k+Q ~ ^fc iQa 

Q-2fc Q + 2A; 



^fc+2Q - c 2iQx 



(33) 



Here, we have used the "pre- averaging" Kubo formula because there is interesting macro- 
scopic spatial dependence in the conductivity. Unlike atomic fluctuations, we do not want 
to get rid of this spatial dependence by averaging over x. For 5 = 1, the above expression 
simplifies further due to the fact (0a) vanishes. Lastly, from the dynamical conductivity, 
we can obtain the optical conductivity by simply setting k = 0. 
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2. A / 1 



For A ^ 1, we can solve Eq. ffTUl) at the high-frequency regime to > Q, [(Oi)^] 1 ^, giving 
the optical conductivity 

,,(„, * = 0, ,) » 1 + tpgt _ MM C08 Qx _ C08 2Qx , (34) 

where 

o; = (2 2(1 - A) r(2A-l)(G A )(°) 2 )^ , (35) 

for 8 ^ 1. Again, this simplifies further for 5 = 1 as the odd modes of the order parameter 
vanish. 

We can also solve Eq. [19] in the low frequency regime to get 

TL 

a y (u, k — 0, x) ~ i— , (36) 
where n s is the superfluid density. Explicitly, 



n s _ x | (0a) (O) 
w 



2 



4(O a )« ( Jl + ^-ljcosgx+(0 A ) (2) Ul + -^f-l )cos2Qx 



(37) 

for 5 7^ 1. As before, for 5=1, the odd modes of the order parameter vanish. 

Using the Kramers-Kronig relations, we deduce the real part of the DC conductivity, 

%t[a y {u = 0,k = 0,x)\ = im e 8{u) . (38) 

We note that the superfluid density vanishes when the scalar order parameter vanishes 
(because ojq = 0), which corresponds to the onset of superconductivity. 



B. Conductivity along the direction of the stripes 
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As noted in Ref . [9] , because the inhomogeneity is in the x-direction, applying an electric 
field in the re-direction sources other independent perturbations even at linear order, making 

9 



the computation of a x complicated. The perturbations that need to be turned on are 

cLud 2 k 



A x - 
A t -> A t 

X = 



(2vr)= 



da; d 2 /c 
did d 2 k 



a t (z;u,k,x) e ^ kx - ut \ 



duj d 2 k 



f 2 (z;u,k,x)e« kx -^. 



(39) 



Here, adding the small perturbation in term of the bulk vector potential sources perturba- 
tions for the bulk electrostatic potential and for both the real and imaginary part of the 
bulk scalar field. We have also let the Fourier coefficients a x , at, j\ and f 2 to retain resid- 
ual x-dependence due to the fact that the subleading solutions are related to the leading 
order solutions whose momenta differ by multiples of Q. This is exactly what we saw in the 
previous subsection (see Eq. 125]) which ultimately leads to the spatial dependence of the 
conductivity. 

The equations of motion are then given by 



(-Wi) 



~h 



[ik + d x ) 



fi 



A 2 t , A(3-A) 



h 



z' 



fi + T A t (Va t + ujf 2 ) = 0, 



z 2 d z - d z f, 



^- + (ik+d x ) 2 



A 2 , A(3-A) 



h 



2iu 
~h 



iUJ 



-j^ a t + (ik + d x ) a a 



2d x ty a x = 0, 



hd 2 z + (ik + d x 



a t -^ 2 a t - 2* A t fx ~ uj * / 2 + icoz 2 (ik + d x ) a x = 0, 



f 2 d z V-Vd z f 2 + z 2 



IUJ 



d z (ik + d x ) a x + — d z a t 
h 



z 2 d z (hd^ 



cj 2 z 2 



IUJZ 



a x -V 2 a x - d x V f 2 + ^> (ik + d x ) f 2 — (ik + d x ) a t = 0. 



h 



(40) 



Even though it is difficult to proceed for a general value of A, we are able to make progress 
in the case of A = 1. The system has a solution which is pure gauge, given by 



at 



fi = 



and 



/ 2 = *. 



(41) 



At large Q, we can neglect A t and expand around z = so that h « 1 and fa 1. We 
then have f\ = and as before, we can evaluate the solution up to the sub-leading n — 1 
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and n = 2 modes 



,(0) 



aW e -iQx 



a t = a\' +a\' e v + a 



( /f + if e iQ ^ + if e"** + . 



+ a (2) e i30x + s (2) g-fiQ* 
(2) o i2Qx _|_ -(2) e -i2Qa; 



F(l) „-»Qa 



(2) o i2Q:z 



+ / 2 (2) e - i2Q:c 



(42) 



At zeroth order, we have the system of equations 



(o) 



iuaf ] + ika® = 0, 



(Oi) 



(oy 



n 



uka^ = 0, 



,(°) 



which is solved by 



At next order, we have 



,(o) 



,(0) 



CO 



2 Al) 



dif 2 



2M) 



d z a t 



(O^ + ik + Qf 



to 2 - (k + Qf ] / 2 (1) - i{Ot)^ [ua? + (k + Q) a« 
-(d)^ ^/ 2 (1) + i (A; + Q) d z aW + 

wr«'<Oi> (0) / 2 (1) + (* + 0)flf > ' 



+ k - (^i) (0)2 l 4 1} + (k + o) k^) (0) /f + 



(43) 
(44) 

^(0i) (1) 4 O) , 
0. 

(O 1 )(°)(O 1 ) (1) «l 0) 
(Oi>W<Oi>«af), 
(45) 



with boundary conditions /. 



(i) 



at z = 0, while at the horizon, and g4 



have the form of an incoming wave, and vanishes. The second equation then becomes 

-<Oi)(°> /f + z (A; + Q) a« + ^a| 1} = . (46) 



Using this, we see that the remaining three equations decouple, 

L k+QJ2 



n 2 ' m 

Jit 



rv 



2„W 



(1) _ 



a 2 a x +ii fe+ ga x — e 



(47) 



11 



The solution that satisfies the correct boundary conditions is 

,<!)_ *(<?1> (1) 



k(Q + 2k) 



) 



uQ(Q + 2k) \ sinfi fc +Q 

(i) _ (Ci) (0) (d) (1) / iQ fe+QZ _ jn k z\ 
x " kQ(Q + 2k) 16 >■ 



sin f2 



(48) 



The barred n = 1 mode functions are obtained from the above expression by replacing 
Q — > —Q, while the n = 2 mode functions and barred functions are obtained by replacing 
->• (Ci)( 2 ) and Q -)• ±2Q. 
Lastly, we need to combine the above solution with the pure gauge solution to ensure 



that for the full solution, at vanishes at the horizon 
solution becomes 



13j . In particular, the leading order 



/i = 0, 



i(Oi) 



(o) 



UJ" 



(o) _ g - e 



Co' 



a (o) = 1 ( e ia kZ _^ e iu k 



k \ uj* 

The conductivity along the direction of the stripes is then given by 



(49) 



a x (u,k,x) 



UJ 



UJ' 



— k 2 e iClk 



Q 



Q-2k 



Q + 2A; 



(Cl) (0) (Cl) (2) ^ ^fc-2Q - ^fc -2iQx , ^fc+2Q - ^fc 2iQx 



4Q 



Q-k 



+ Q + k 



(50) 



We note that at low temperature, for A = 1, the optical conductivity is isotropic. 



IV. NEAR THE CRITICAL TEMPERATURE 



As we have explained previously, the analytical techniques developed in the previous 
section are valid when the order parameter is a lot larger than the temperature. In order to 
study the system near the critical temperature, where the scalar order parameter is small, 
we need to be more careful in handling the boundary condition at the horizon. We can make 
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progress in our analytical study by considering the long wavelength, small frequency regime 
k,u < 1. 



A. Conductivity perpendicular to the direction of the stripes 

For a y , for convenience, let us factorize the appropriate boundary behavior 

Ay = J ^^^k)(l-zr^e^\ (51) 

so that a y is regular at the horizon. Furthermore, we express the perturbation as a Taylor 
series in u ~ k 2 ~ £ 2 <C 1 

a y = a y o + a yu; u + a yk k 2 + a yi £ 2 + . . . (52) 

where £~ (O A ) (0) - 

Moreover, each field we divide into the n = mode and the n — 1,2 modes that are 
suppressed at large Q (higher modes are even more suppressed). As before, we shall use a 
superscript to denote this hierarchy. 

At leading order in Q, we have 

(1 - z 3 ) d 2 a° y0 + ^d z a° y0 = 0, (53) 
whose only solution that is regular at the horizon is the constant 

= = 0) = a° y (z = 0). (54) 

At O(oj), we then have 

o ( s i o, m r dz ' r 1 ^ + z" - 2(z")\ „ .... 

= 3 ^ = 0) J W) I 1 - z» dZ ' (55) 

while at 0(k 2 ), we have 

a° yk (z) = -al(z = 0)^^ r) j\z". (56) 

Finally, at 0(£ 2 ), we have 

,0 („\ _ „Qf„_ n \ / az / AJ'J' 2 ^- 2 



<&(*) = -a y (z = 0) ^ / dz" V' 2A " 2 [F<V)] • (57) 
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At sub- leading order in Q, all of the fields vanish except 
dz' rl 



alz(z; k) 



K Z ') Jz> 

dz' rl 



o Kz') 



dz „ z „2A-2 F{0)F(1) 



dz « z »^-* F (0) F (2) 



a y (z — 0;k — Q) + a y (z = 0; k + Q) 
a y (z = 0;k- 2Q) + a y (z = 0;k + 2Q) 



(58) 



The dynamical conductivity perpendicular to the stripes is then given by 

2 



k,x) = l +— + L ^ A; J / dz-^ x \F^ + 2F^cosQx + 2F^cos2Qx] , 



(59) 



from which we deduce the real part of the DC conductivity fl38|) with the superfluid density 
given by 

•i i?(o) 



n 



s..V 



[(C A ) (0) ] / ^4=2A [ F(0) + 2F(1) COS ^ X + ^ COS 2 ^ X ] • 



(60) 



B. Conductivity along the direction of the stripes 



To obtain the conductivity along the stripes, a x , the perturbations that need to be turned 
on are 

dw d 2 k 



A 



a x (z;u,k,x)(l- z)~ iul/3 e 1 ^"^ , 



(2-Kfk 



* ->■ * + 



(2tt) 3 



* (0) (2) /a (z; u,k,x) (1 - 2)-^ /3 e i(to - w *) , 



X 



(2tt) 3 

^ /a(z;Wj kjX) (1 _ ,)-/3 e .("). 



(61) 



Again, we have conveniently taken out a factor which captures the behavior at the horizon, 
so that all functions should be regular at the horizon. 
The equations of motion for the n = modes are 



2iu 



[l + z + z 2 



2 + z 3 



3 

iui(2 + z) 2iu 



2h 



(o) 



fl + Z + ^ (0) ^ + 2 )( 4 + ^ + ^) ,2 



9(1 + z + z 2 ) 
--~A t (af + iojf. 



n 



(0) 



(or 

2 y > 



(62) 
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liuj . o, 2 + z 3 
— (1 + z + z 2 ) 



2h 



\£(o) 



iu{2 + z) 
3z 



2iu 
~3 



^(l + ^ + ^ 2 



2iw 



9(1 



(0) 



(63) 



2iw 



k 2 + 



2(^W) 2 w (3i-w)(l + z + z ) 



2(^(o)) : 



9(1 -z) 



2A*f® + iw/f 



,(°) 



(o) 



(64) 



3(1 -z) 
2iu> 



,(0) 



2(^(o)) : 



3(1 

(l + 2 + ^ 2 )-32 2 <9 2 af 



iw, N 2 8 + 6^ + 3^ 2 + ^ 3 2(tf(°)) 

y (1+2 ^ )+ " 9(1 + , + , 2 ) V 

uk 2 



d f (0) + - 

+ 3(1-^ 



f (0) 



(65) 



,(0) 



x 

(0) 



h 



2^(^Yff ) . (66) 



Here, we have neglected the n — 1 mode of A t , because it is exponentially suppressed at 
large Q. Thus, 

A « /i(l - «5)^ 0) « M (l -5)(l-z)= Ml " 2) • (67) 

Since u, k and \& are small, we can express the modes as Taylor series in w, /c 2 and £ 2 as 
before, 



/r 


_ f(0) 
~~ J 10 




+ A ( fe 2 


f (l) 

J2 


_ f(0) 
— J 20 




+ / 2 ( fe 2 




(0) 
— a i0 


+ aS° ) w + aJ° ) /c 2 


+ < } e 2 


<£> 


-a (0) 

_ a x0 


+ 4°> + ai^ 2 





(68) 

We shall solve the field equations starting from the horizon, so that all integration constants 
will be of zeroth order there. The boundary values will then be determined in terms of these 
horizon parameters. 
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At zeroth order, we obtain from Eqs. (164]) . (1651) and (1661) . 

#*2? = ^offi = 0, (69) 

which is solved by 

a to( z ) =oto(l-z), , a£° (*) = a . (70) 
The two parameters are related by demanding regularity of ai2 in Eq. (1631) 

a = -ia , (71) 



and thus, we obtain 



From Eq. (1661) . we also obtain 

«2W = 0, , aS ) (,)=-o£^£^V^- 2 [^)(^)] 2 . (73) 
From Eq. (1651) at 0(£ 2 ), or Eq. (I63p at zeroth order, we then deduce 

f°»(z) -B a r i:J Z ' 2 ~ 2A f 1 dz" (74) 

where the integration constant /?2o will be determined later. 
Next, let us take another look at Eq. ( 1641) to obtain 

a%(z) = --ff M d *'i dz \ + y +z ,* > a< $^ = «°/ dz 'J zl 1+ d J +z »> - ™ 



Now, we can solve Eq. (162|) at zeroth order, 

= p 10 + 2« f dz' , /' dz" MZ,,) 2 ^yj} , (76) 

10 KJ Jo h(z')[FM(z')} 2 L> l + z" + z" 2 z" 2 - 2A 1 ; 

where the integration constant j3 w will also be determined shortly. We can then obtain 



Z n2-2A \ + z " + z " 

from Eq. ( 1641) . 

At 0(u), Eqs. (J62J) and fl63]) yield 

M' )= ^ + L"mm®Lm^' (78) 
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for i = 1,2, where 



F 1 (z) = 2A t (a$+if. 



(o) 

20 



+ -h 

3 



2(1 + z + z 2 )d z f$ 



F 2 (z) 



-i a 



'70 



+ ha<g>) - 2tAJ 



(0) 
10 



2(i + ^ + ^ 2 )a 2 / : 



(o) 

20 



2 + 2 



2(1 + 2 + 2 



2\ ^2 



(0) 
20 



(79) 



We are now ready to determine the integration constants. This is done by forbidding the 
bulk scalar field to generate source terms for the scalar order parameter of the boundary field 
theory. In other words, we are going to determine the integration constants by excluding 
the unwanted behavior / 1)2 ~ Z 3 ~ 2A at the boundary (2 = 0). Here, we need to distinguish 
between two cases, 1 and A = 1. 

For A 7^ 1, to isolate the unwanted behavior, we write 



(80) 



where 



Ao = /?io - 2a / dz' 



/2-2A 



h(z>) [FW(z')Y Jo 



f w = 2a / dz' 



/2-2A 



*' „ A(z") [F^(z")Y 

aZ l + z » + z «2 2//2-2A > 
2 



1 „ A t {z") [F«»{z"j\ 

y" _L „//2-2A 



(81) 



'0 h(z')[FW(z>)} z Jo l + z" + z" z z" 

Hereupon, we use the notation * and " to describe the wanted and unwanted near boundary 
behavior, respectively, so that g ~ z° while g ~ z 3 ~ 2A . 
Similarly, for / 2 q, we have 



fw — Ao + /: 



20: 



where 



Ao = #20 + "o / dz' 



z 



/2-2A 



h(z>)[F(°)(z>)Y Jo 



* [F<«V)]" 



z 



,2-2A ' 



Ao = -"o / dz' 



/2-2A 



Next, we split 



h(z')[F(°)(z')Y Jo 



z 



//2-2A 



(82) 



(83) 



(84) 
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for i — 1,2, with 



2(1 + 2 + + ( - 2±£ + 2(1 + z + z 2 ) 9 -^ ) 



2(1 + 2 + 



(o) 

10 



2 + 2 



+ 2(1 + 2 + 2 



$(0) 



f(0) 
J 10 



and 



7i = -i (aff + fcaffi) - 2iAj. 

2(1 + z + z 2 )dj§ 



(0) 
10 



3 



l±i + 2(1 + 2 + ^ )F . 



20 



= -2iA t f$ + i/i 



2(l + 2 + 2 2 )a 2 / ; (0) 



20 



2 + 2 



2(1 + z + z 2 )-^j— J f. 



(o) 

20 



Therefore, 



/iw(z) = / dz' 



z' 



/2-2A 



1 



z 



/2-2A 



(85) 



(86) 



1 AJI 



dz" Tiiz"\ 



o ~~ h( Z i)[F(°)(z')} 2 J Z ' Hz") z" 2 ~ 2A J h(z>)[FW(z')f Jo h(z") z" 2 " 2A ' 

(87) 



At lowest order, we then obtain 



2iCan ,„ iVan 

A(3 10 + BB 20 = , A8 20 -B8 W = 



to 



(88) 



where 



2 + z 



+ 2(1 + 2 + 2 2 



A=- dzz 2A ~ 21 
3 Jo 

dz- % \F®] 

1 + z + z 2 1 J 



Therefore, 



where 



8 

00 



-1 z 2A ~ 2 
B = 2fi 8 I dz - ■ , 9 , (89) 



2 

1 + z + z' 2 



V 



dzz 2A ~ 2 [F^] 2 . (90) 



(91) 



2AC + BV 



x 2 



AV - 2BC 



(92) 



A 2 + B 2 ' z A 2 + B 2 ' 
For A = 1, /3jo (i = 1,2) are found by demanding the absence of a z 2 term in the near 
boundary asymptotics of ^ and This implies 



^/l,2 + -3-/1,2 = 



(93) 
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at z = 0. We then have 

r 



A 



W + Bf3 20 



2iCa 



A 



1 



/3 20 - B/3 



War 



10 



where ^4, £>, C and X> are evaluated at A = 1. Therefore, 

2 [A + |) C + BV {A+\)V- 2BC 



Xi 



%2 



(94) 



(95) 



At a glance, it seems that there is a discontinuity at A = 1. However, this is not true. To 
see that we do in fact obtain the correct limit, let A = 1 + e. Then, we have 



A 



3 Jo 



dz z 



2c 



2 + z 



z 



2(1 + z + z 2 



1 + e d z F® 



z 



3 Jo 



dz z 



2c 



2e + z + 2z 2 . 2 .d z F 
- + 2(1 + z + 2 2 )— - 



0(6) 



3 Jo 



dz 2ez 2e ~ 1 + .4 
(9(e) 



A=l 



0(e) 



A=l 



(96) 



showing that we recover the expression used in the A = 1 case as e — > 0. 

Next, let us consider the n — 1 modes. For convenience, rescale -)• /^/A:, / 2 (1) 
frp /k and — >■ dp /k. Then the equations of motion are 



2 Ai) 



+ 



2zuk 2N 2 + z 3 o ,<9 2 # (0) 



d z f\ 



(i) 



ia;(2 + ^) 2^ frgP) ^ + 2)(4 + z + , 2 ) 

3^ + 2 + 9(1 + z + z*) (* + C ^ 



/J 



(i) 



2iu 



(1 + z + z 2 ) 



2 + z 3 



2h 



h 



,(1) 



t u i 



2F(°) a * 



(0) 



iuf. 



(i) 



(97) 



+ 



33 



9(1 + z + z 2 ) 



f: 



(i) 



2iuj 



A t f[ l) +i\^-af ) + (k + Q) a« 



7 Co' 



+ ^ 



F 1 , ^ 

iO=—a {0) 



hd 2 z a^ + 2 -^(l + z + z 2 )d z a^ 



(k + Q) 2 + 
4fc^(°)^( 1 ) r (0) 



2 

2(*(°)) z w (3i-w)(l + z + z ) 



(98) 



9(1 -z) 



(o) 



(i) 



(o) 



(99) 
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i(k + Q) d e a,W - 



uk 



,(i) , ™n 



3(1 -z) 
2(^(0))' 



-d z a\ 



3(l-*)/T* 
/ 2 ( 



_JW_,(1) 

3(1 -z) i2 



+ 



2$(°)^( 1 ) 



0,/ 2 (1) 



3(1 - z y 2 



/ 2 (0) <9^ ), (100) 



— (1 + 2 + 2 2 ) -32 2 
o 



+ 



^ 8 + 6^+3^^ 2(tt(°>)- 

y (1+22)+u; 9(1 + ,+^) pr- 

_w(fc + Q) (1) 2i(fe + Q) 



(^(o)) 2 /: 



(i) 



As with the n = modes, we expand the n = 1 modes as 



(o) 

x ■ 



(101) 





— J 10 


+ /£>«- 




+ A?? 2 


f(i) 

J2 


_ f (i) 


+ - 






a? 


-a (1) 
— "to 








<£> 


-a (1) 


+ + «2* 2 





(102) 

in the small oj, k and £ regime. Then Eq. (199]) . along with the boundary conditions 0^(2: = 
0) = and a$(z = 1) = 0, give us a$(z) = 0. Similarly, af](z = 0) = a^}(z = 0) = 0. Eq. 
f llOOl) . along with the boundary condition ){z = 0) = 0, and demanding regularity at the 
horizon, gives us a^(z) = 0. Similarly, a£!(z = 0) = a^(z = 0) = 0. 

In order to solve Eq. ( |98|) at leading order, we do a coordinate transformation z = Qz, 
and keep the leading terms in Q to obtain 



6-,$) -f. 



(i) 

20 



(0a) (1) 



O: 



(103) 



The solutions to the homogeneous part of this differential equation can be expressed as linear 
combinations of Bessel functions 



- 3-2A , 

Z 2 J 2A-3 \ —lZ) 



. . 3-2A 

and z 2 I2A-3 (—iz) 



(104) 
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However, the latter will give rise to a source term at the boundary and thus, the solution 
that satisfies all necessary boundary conditions in the original radial coordinate is given by 

a (Oa) (1) 



f (i) _ 

720 Q (O a )C°) 



M-iQz) 



where 



Ja(u) 



r A 



3-2A 

u\ —5— 



From Eq. fllOip . we then obtain 

dz' rl 



(i) 



o Kz') 



dz" z" 



( F (o)y 



J2A-3 (it) 



F« (Ca) (1) 



(105) 



(106) 



(107) 



At zeroth order, Eq. (1971) . along with the boundary conditions, and demanding regularity 
at the horizon, yields = 0. Substituting this into Eq. ( )99|) at 0(£ 2 ), we get af^ = 0. 

Moving on to the other modes, the n — — 1 modes are similar, and the n = ±2 modes 
are obtained by substituting Q 2Q, (C A ) (1) -> (0 A ) (2) and -> F^ 2 ) into the results 
for n — 1 modes. 

The boundary values of the perturbations of the gauge fields are then given by 



at(0) « a + 0^(0)0; + a^'^k 2 + a^(0)£ 2 « a 



1 + — «0 A >(°>) : 



a s (0)«a + 4°2(0)u; + 4>)£ 2 



where 



(108) 



(109) 



ft z' 2 " 2A 

Clearly, this is not a general solution, because the boundary values of a t and a x are not 
independent of each other. To correct this, we need to mix in an arbitrary amount of pure 
gauge solution 



a t (l - z 

We then have 

a t (0) = a 



-ito/3 



-iu^o and a 



1 + ^((O a )(°)) ; 
ui 



; (l - Z )-*» /3 = ^ 2 7o . 



flx(0) = -z'a + ^fc 2 7o- 



(110) 



111) 



Solving for the arbitrary parameters, we obtain 
k 2 a t (0) +ua x (0) 



k 2 



1 + 



2iC 



((Ca) (o) ; 



To 



a t (0) - 


- z 


\ + MZ«0 A )(o)) 2 ~ 


a*(0) 


fc 2 


y + MZ((O A )(0)f 







(112) 
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We can now read the retarded Green functions from the on-shell action. The (charge) 
density- density retarded correlation function is given by 



[pp)(u,k,x) 



k 2 (cu + 2tC[(0 A }(°)] 2 ) -iuj 



7T 



UJ 



2ip s x 1 [(0 A )W] 2 [ - 
Jo 1 



2 f 1 dz [F^\ 



z + z l z 



2 ~2-2A 



with poles at 



UJ 



1 ± \ 1 -AC 



k 2 



When (Oa)^ <C k, we get a diffusion pole 



, (113) 



114) 



uj = —ik^ 



(115) 



showing that the behavior of the system is closer to a normal fluid then a superfluid, while 
when (Oa) ^ 3> k, the behavior is closer to that of a superfluid and we get poles at 



UJ 



±Vc{o A )®k, 



(116) 



which is identified as the second sound with the speed of the second sound given by 



v = Vc (O A } {0) . 



;ii7) 



We can also obtain the thermodynamic susceptibility, 



d 2 P = \im(pp)(0,k, 



ifc->0 



X) 



Psx, f 1 dz [f(°)] 2 
C Jo 1 + z + z 2 z 2 ~ 2A 



(118) 



which is positive definite. 

Next, we can read the current-current retarded correlation function from the on-shell 
action, and then take the limit k — > to get the optical conductivity in the direction of the 
stripes, 



a x (oj, k = 0, x) = 1 + i—^- + . . . 

UJ 



(119) 
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where 



r 1 p(o) 

[(o A ) (0) ] / dz ~i^K i F{0) + 2p[l) cos Q x + 2p{2) c °s2gx] 

Jo z 

+ 4(O a ) (0) (Oa) (1) / dz^^M-iQz) cos Qx 
Jo z 

+ 4(0a) (o) (Oa> (2) f 1 d z l^-^j A (-2iQz) cos 2Qx . (120) 
Jo z 



It is easy to check that n s is real. It is related to the DC conductivity via 

Comparing with the conductivity perpendicular to the stripes ( 1601) , we see that the 
optical conductivity is anisotropic, 

<j x (uj, x) — a y {co, x) = i — ■ — + . . . (121) 

where 

2 



4(0 A > (O) (O A > (1) f dz^^J A {-iQz)cosQx 
Jo z 

+ 4(O a ) (0) (Oa) (2) / dz^^J A (-2iQz) cos 2Qx . (122) 
Jo z 



V. CONCLUSION 



In this article, we have calculated the conductivity of a strongly coupled striped super- 
conductor at large modulation regime, both in the directions perpendicular and parallel to 
the stripes. In the dual gravity picture, this is done by adding perturbations in term of 



io| 



vector potential on top the solution of Einstein-Maxwell-scalar found in Ref. 

We found that below the critical temperature, the imaginary part of the optical conduc- 
tivity has a pole at u — 0, which implies that the real part of the conductivity contains a 
term that is proportional to 5(u), thus establishing that this model indeed exhibits super- 
conductivity below the critical temperature. 

At T < T c , we found that the conductivity is anisotropic and it exhibits 1/u behavior 
at small frequency. At low temperature T <C T c , we were only able to calculate a x ana- 
lytically for A = 1, and in this case we found that the optical conductivity is isotropic. 



As the anisotropy at low temperature has been measured [l4j, our results seem to rule out 



A = 1. It will be interesting to study this model at low temperature for different values 

23 



of A numerically and compare it to the behavior observed experimentally at large 

frequency. 

It will also be interesting to study the anisotropy at T > T c by introducing fermion 
into this holographic set-up and see whether there is something akin to the "dimensional 



crossover" as studied in 
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